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Abstract 

Recently the AMS-02 collaboration has published the measurement of the cosmic antiproton to 
proton ratio p/p and the p flux with a high precision up to ~ 450 GeV. In this work, we perform a 
systematic analysis of the secondary antiproton flux generated by the cosmic ray interaction with 
the interstellar gas. The uncertainty of the prediction originates from the cosmic ray propagation 
process and the hadronic interaction models. Although the cosmic ray propagation parameters have 
been well controlled by the AMS-02 B/C ratio data for a specified model, different propagation 
models can not be discriminated by the B/C data. The p flux is also calculated for several hadronic 
interaction models, which are generally adopted by the cosmic ray community. However, the results 
for different hadronic models do not converge. We find the EPOS LHC model, which seems to fit 
the collider data very well, predicts a slightly lower p/p ratio than the AMS-02 data at the high 
energy end. Finally we derive the constraints on the dark matter annihilation cross section from 


the AMS-02 p/p ratio for different propagation and hadronic interaction models. 


PACS numbers: 96.50.8-,95.35.--d 


I. INTRODUCTION 


Antimatter searches in cosmic rays (CRs) are especially important for astrophysics and 
dark matter (DM) indirect detection. In recent years, there were great progresses in the 
measurements of CR antimatter particles. In particular, the Alpha Magnetic Spectrome- 
ter (AMS-02) launched in 2011 has provided unprecedentedly precise CR. measurements. 
Based on the results from AMS-02 |1] and previous experiments, such as PAMELA |2] and 
Fermi-LAT|3], the properties of CR antimatter particles have been extensively studied in a 


quantitative way in the literature. 


Recently, the AMS-02 collaboration has published the measurements of the p/p ratio and 
the p flux up to ~ 450 GeV with a high precision |4]. The major challenge to interpret the 
AMS-02 data and search for potential DM signals is how to predict the secondary antiproton 
flux precisely. The secondary antiprotons are generated by collision between the high energy 
cosmic rays and the interstellar gas. Its flux is determined by the CR propagation process 
and its hadronic interaction with the interstellar gas, both of which have relatively large 
uncertainties at present. In this work, we perform a systematic analysis of the secondary 
antiproton flux considering these possible uncertainties carefully. Some relevant discussions 


can also be found in Refs. |5H14! 


The CR propagation in the Galaxy involves many complicated effects, such as the dif- 
fusion, energy loss, convection, and reacceleration. The propagation parameters describing 
these effects can be determined by fitting to the secondary-to-primary nuclei ratios, such as 
B/C and (Sc+Ti+V)/Fe, and unstable-to-stable ratios of secondary nuclei, such as !?Be/?Be 
and ?6A]/?' A] 


collision with the interstellar matter when they are propagating. However, the current mea- 


, where the secondary nuclei are generated by the primary nucleus 


surements of these ratios remain not sufficient to distinguish different propagation models. 
Therefore the propagation is one of the main sources of uncertainty when predicting the 
secondary antiproton flux. Another important uncertainty of propagation comes from the 
solar modulation lid. This effect significantly affects low energy CR spectra, but is difficult 
to precisely quantify. 

In this work, the GALPROP package is used to solve the CR propagation equation 
" We adopt three kinds of propagation models, namely the diffusion-reacceleration (DR) 


model, the modified diffusion-reacceleration (DR-2) model, and the diffusion-convection 
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(DC) model. Propagation parameters are determined by fitting the available B/C data 
with a Markov Chain Monte Carlo (MCMC, m algorithm , which is powerful for 


surveying high dimensional parameter space. 


The hadronic interaction between CR particles and interstellar gas is another important 
source of uncertainty for calculating the antiproton flux. Although Quantum Chromody- 
namics (QCD) is well-established as the theory of strong interaction and has been confirmed 
at collider experiments, only the processes with large momentum transfers can be predicted 
by the perturbative calculation from the first principle. It is not possible to calculate the 
forward scattering processes with multiparticle production by QCD, such as the CR inter- 
action with interstellar gas considered here. In order to calculate such hadronic processes, 


some simplified assumptions and phenomenological/empirical parametrization are needed. 


In the cosmic ray study, it has been a long time difficulty to deal with the hadronic 
interaction. Many phenomenological forms have been proposed in the literature to describe 
the hadronic interaction. Unfortunately, there is no consensus on the empirical hadronic 
model so far. In this work, we discuss the impact of hadronic models on the antiproton 
production. We find that the hadronic interaction induces the largest uncertainty in the 
prediction of the antiproton flux. 

Based on the predicted p flux and the AMS-02 measurement, we derive constraints on the 
DM annihilation cross section in different hadronic interaction and CR prorogation models. 

This paper is organized as follows. In Sec. II, we describe the CR propagation processes. 
In Sec. III, we introduce the hadronic interaction models and compare the prediction with 
accelerator data. The secondary antiproton flux and p/p ratio prediction are presented in 
Sec. IV. In Sec. V, we investigate the implications for DM annihilation from the AMS-02 
p/p measurement and the CRs prediction. Finally, we give the summary and discussions in 


Sec. V. 


II. PROPAGATION OF GALACTIC COSMIC RAYS 


After accelerated in sources, Galactic CRs are injected and diffuse in the interstellar space 
and suffer from several propagation effects before arriving at the Earth. A conventional 


assumption is that CRs propagate in a cylindrical halo with a half height z}, beyond which 


CRs escape freely. The propagation equation can be expressed as 
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where Q(x, p) is the CR source term, v = v(x,p,t) is the CR density per momentum 
interval, p = dp/dt is the momentum loss rate, and the time scales 7; and 7, characterize 
fragmentation processes and radioactive decays, respectively. 

It is generally believed that CR particles are accelerated in supernova remnants (SNRs). 


Thus the spatial distribution of the CR sources follows the SNR distribution as 


foo) o 


where rg = 8.5 kpc is the distance from the Sun to the Galactic Center, and z, ~ 0.2 kpc 
is the characteristic height of the Galactic disk. We adopt a = 1.25 and b = 3.56 in 
the calculation. The injection spectra are assumed to be broken power-law functions with 
respect to the rigidity. 

The spatial diffusion coefficient Dys can be parametrized as 


Dzs = Dop” (R/ Roy, , (3) 


where R = pc/Ze is the rigidity, 3 is the CR particle velocity in units of the light speed c, Ro 
is the reference rigidity and Do is a normalization parameter. Although the slop of diffusion 
coefficient ô is predicted to be 6 = 1/3 for a Kolmogorov spectrum of interstellar turbulence, 
or ó = 1/2 for a Kraichnan cascade, this parameter is usually treated as a free parameter 
when explaining data. The factor of 3” denotes the effect that the diffusion coefficient could 
be altered at low velocities due to the turbulence dissipation. We take 7 = 1 as the standard 
diffusion plus reacceleration (DR) case taken in GALPROP. We refer the case with 7 = —0.4 
as the DR-2 model, which can improve the secondary fitting as discussed in [29]. 

The CR reacceleration process due to collisions on interstellar random weak hydrody- 
namic waves can be described by the diffusion in momentum space with a coefficient Dpp, 
which is related with D,, by 


Ap^v? 
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where v4 is the Alfvén velocity and w is the ratio of the magnetohydrodynamic wave energy 
density to the magnetic field energy density. Since Dpp is proportional to vÀ /w, we can set 
w = 1 and just use v4 to characterize the reacceleration effect. The convection velocity V, 
is usually assumed to linearly depend on the distance away from the Galactic disk and the 
convection effect can be described by the quantity dV./dz. Therefore, the major propagation 
parameters involve Do, 6, Ro, n, va, dV./dz, and zp. 

When CRs propagate in the solar system, their spectra with R < 20 GV are significantly 
affected by solar winds. This is the solar modulation effect, which depends on solar activities 
and varies with the solar cycle. It can be described by the force field approximation 
with only one parameter, i.e., the solar modulation potential $. 

Due to fragmentation and radioactive decays, secondary CR particles are produced in 
propagation processes. As a result, secondary-to-primary ratios and unstable-to-stable sec- 
ondary ratios, such as B/C and "Be/?Be, are sensitive to propagation parameters, but 
almost independent of primary injection spectra. Thus the measurements of these ratios are 
useful for determining the propagation parameters (see e.g. Refs. bi j 

Note that there are degeneracies between the propagation models with the reacceleration 
process and those with the convection effect. Current measurements are not sufficient to 
distinguish these models. Therefore, we separately consider the DR, DR-2, and DC models 
when calculating the antiproton flux. Ro is taken to be 4 GV in the DR and DR-2 models. 


In the DC model, however, Ro is set to be a free parameter and ô = 0 is imposed for R < Ro. 


III. HADRONIC INTERACTION MODELS 


Although QCD has been proved to be a very successful theory of the strong interaction, 
we do not have a realistic technique to calculate the bulk of multiparticle production in 
hadronic interactions from the first principle. Many phenomenological and empirical models 
have been developed to explain the collider and comic ray data. These models include a 
large number of parameters that encode the fundamental physics and phenomenological 
descriptions of the fragmentation process. 

These hadronic models could be classified into three categories according to their pur- 
poses. HERWIG 33) PYTHIA[33 and SHERPA 4 are generally adopted in high energy 


physics (HEP) studies, and focus on hard-scattering processes at accelerator experiments. 
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In contrast, QasJETO1[3. QGSJET mM and SIBYLL are designed for simu- 
lating the extensive air showers caused by the high energy CR particles. Therefore, these 
models focus on the bulk production of soft particles with low transverse momenta, and 


emphasize on providing a reasonable extrapolation for higher energy and wider phase-space 
ld fall in between these two cate- 


regions. PHOJET|4 , DPMJET|M and Eposs 
gories. These models can also approach the HEP models regarding the hard process, and 
are adjusted to well reproduce the low energy accelerator data. 

The HEP models listed above are not considered in our calculation, as they are not 
suitable for dealing with CRs interaction with the interstellar gas. As the hadronic models 
proposed for the CR interaction mainly focus on the high energy region, the description of 
the interaction for CRs below 100 GeV may be not precise. Fortunately, in recent years many 
new data sets measured at LHC and fixed-target experiments are available. New versions 
of QGSJET II-04 and EPOS LHCl4G have been released whose model parameters 
are tuned to fit these data » . The EPOS LHC model fits the available accelerator 
antiproton data in all energy regions. Recently the QGSJET I-04 has been slightly modified 
to better reproduce the accelerator antiproton data at low energies, referred to as QGSJET 
II-04m here . In the following, we use EPOS LHC and QGSJET II-04m to calculate the 
CR antiproton flux. Another version of EPOS, namely EPOS 1 99 [4] is also considered 
in the calculation, since this version has been more carefully tuned to fit the low energy 
accelerator data, such as the measurements of NA49 and NA61 [55]. 

For comparison, we also consider the default hadronic model embedded in GALPROP 
, which is referred to as Tan & Ng + BPO1. In this model, a parametrization given by 
Tan & Ng [57] and a code from Barashenkov & Polanski are combined to estimate the p 
spectrum. Other models for high energy air showers (QGSJETO1, QGSJETII and SIBYLL) 
are not used to predict the CR antiproton flux, as they do not aim at describing the low 
energy CR interactions considered here. 

In the following we compare the expectations of different hadronic models with the ac- 
celerator measurements of the process pp —> pX. In Figs. [I] and 2] we show the cross section 
of pp > p+ X as a function of rp = 2p,/,/s in the center-of-mass system for specified 
transverse momentum pr and integrated cross section over pr for different hadronic models 
respectively. The expectation is compared with the NA49 data Ed 9. which is a fixed 


target experiment with the incoming beam momentum of 158GeV /c. In Fig. [| we consider 
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FIG. 1. The differential cross section of E d?^o /d?p for the p + p — p+ X process compared with 
the NA49 measurement [58], with variables given in the center-of-mass system. Three groups of 
curves from the top to bottom correspond to the pr ranges of [0.05 — 0.15|GeV, [0.45 — 0.55]GeV 


and [0.85 —1]GeV, with a rescale factor of 3, 1 and 0.3, respectively. xp is defined as xp = 2p;/ s. 


1/7 x, do/dx, (mb) 


[7|---- EPOS 1.99 — EPOS LHC 
[E [----QGSJETI-0M4 — QGSJETII-04m — 
| |---- SIBYLL 2.3 


7 larsa drnrnni larraidh 
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 


o 


FIG. 2. The differential cross section of 1/z rgdo/dxp for the p+ p — p-- X. rg is defined as 


Lp = 2Ep/vs in the center-of-mass system. 
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FIG. 3. The differential cross section of E d?c/d?p for the p + p > p+ X process compared with 


the CERN ISR experiment|59] with v/s = 53GeV. 
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FIG. 4. The antiproton spectrum 1/Nd?N/dydp for the p +p — p+ X process as a function of 
the transverse momentum pr in the center-of-mass system with ,/s = 900GeV. The d y is 


cut with the condition |y| < 0.5. Also shown are the data from the ALICE experiment [51]. 


three representative ranges of pr standing for the cases of low, medium and high pr. As 
can be seen that all the models reproduce the cross section integrated over pr well, but only 
QGSJET II-04m, EPOS LHC and SIBYLL could fit the differential cross section for the 
specified pr. 

For a higher energy region, we compare the model prediction with the data from the 
CERN ISR experiment with the center-of-mass energy of \/s = 53GeV. Such energy is 
relevant to the secondary antiprotons with energies of hundreds of GeV. In Fig. [3] we show 


the expectations of differential cross section for pp — pX and compare with the ISR data. 


As can be seen that all the selected models reasonably reproduce the data in low zp region, 
while slightly overestimate the antiproton in the high zp region. 

Finally, we compare the expectations with the measurement of the ALICE experiment 
at ys = ad 
selected hadronic models could give a good fit to the ALICE results except EPOS 1.99. 
Thus we conclude that the EPOS LHC and QGSJET II-04m models are suitable to predict 


. The antiprotons are for the rapidity |y| « 0.5. We find that all the 


the CR antiproton as they can interpret all the accelerator antiproton data well. Note that 
the center-of-mass energy 900GeV corresponding to a CR proton energy of 430TeV is far 
beyond the energy region relevant to the AMS-02 antiproton data. Therefore the EPOS 
1.99 model is also adopted to calculate the CR. antiproton flux, as it can interpret the low 


energy accelerator data well. 


IV. ASTROPHYSICAL PREDICTION FOR THE p/p RATIO 


In this section, we calculate the secondary antiproton flux generated by CRs when they 
propagate in the Galaxy, and study the uncertainties from the propagation and hadronic 
interaction models. 

We use the numerical tool GALPROP to solve the propagation equation. The 
propagation parameters are determined by fitting the B/C and !?Be/?Be data. In order 
to improve the fitting efficiency, the MCMC method is employed to derive the posterior 


probability distributions of propagation parameters. By fitting the B/C ratio data from 
AMS-02 A and ACE ld, and the '°Be/*Be ratio data lal] we obtain the mean values 
and lo errors of propagation parameters for the three propagation models, as shown in 
Table [I| Here the proton injection spectrum is determined by fitting the AMS-02 and 
CREAM proton data. The proton spectrum hardening at ~ 200 GeV can enhance the 
p/p ratio at the high energy end. 

In Fig. [5] we illustrate the p/p ratio in the DR propagation model for different hadronic 
interaction models. The QGSJET II-04 overestimates the low energy antiprotons. This is 
not surprising since it does not reproduce the low energy accelerator data well as shown 
in Fig. The fitting is improved in the modified QGSJET II-04 n" The SYBILL 
underestimates the low energy antiprotons greatly. The reason is that SYBILLL does not 


extend to the cases with the incoming particles energy below ~ 60GeV. We do not consider 
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TABLE I. Mean values and lo uncertainties of the propagation parameters derived through fitting 


the data of B/C and !?Be/?Be ratios in three propagation models. 


DR DR-2 DC 
Dy (1075 cm? s7!) 6.58+1.27 . 3.59 0.88 1.95 — 0.50 
ô 0.333 + 0.011 0.423 + 0.017 0.510 0.034 
Ro (GV) 4 4 4.71 + 0.80 
va (km s^!) 37.8 + 2.7 22.6 + 3.1 / 
dV,/dz (km s^! kpc ^!) / / 4.2 + 3.2 
zn (kpc) 4.71.0 3.5 + 0.8 2.5 -E 0,7 
9pjc (MV) 326 + 36 334 + 37 182 + 25 
n p/p 
x10 
Q 
IQ - 
- + AMS-02 
0.254 | 
0.2 
0.15— 
0.1— 
T d — Tan&Ng + BPO1 "=" SYBILL 2.3 
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E eu sU [e QGSJETII-04 — — QGSJETII-04m 
1 10 10? R(GV) 


FIG. 5. The p/p ratio expected from different hadronic models in the DR propagation model, 


in comparison with the AMS-02 measurement |4]. The propagation parameters used here give 


the best-fit to the B/C ratio, adopted from Ref. [25]: Do = 6.58 x 107? cm? s^, 6 = 0.33, v4 = 


37.8kms-!, zp = 4.7 kpc, vı = 1.81, v = 2.40. The solar modulation potential is set to 800 MV. 


the two interaction models later. 


However, even the other four hadronic models provide diverse predictions. Note that the 


propagation parameters have been adjusted to fit the B/C data. Therefore the discrepancies 
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between the p/p predictions and data may be due to the uncertainties in the nuclear reaction 
cross sections for C, O — B or in the hadronic interactions for pp + pX. The uncertainties 
may also come from the heavy element interactions pA(AA) — pX. In order to include 
such effects, we introduce a free factor cp to slightly rescale the antiproton flux in order to 
fit the AMS-02 data. Therefore we present the prediction of p flux and p/p ratio, which is 
"calibrated" by the AMS-02 data, in the following, since the AMS-02 collaboration provides 


very precise measurements and well-controlled systematic errors. 
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FIG. 6. The p/p ratio (colored) and p flux (gray) for different hadronic models in the DR model, 
comparing with the AMS-02 measurement |4]. The bands indicate the propagation uncertainty in 
9596 C.L., while the line is the best-fit case to the p/p ratio. We have tuned the p/p ratio with 


suitable rescaling and solar modulations to fit the data, and the p flux prediction corresponds to 


this fit. 


In Figs. [6] [7] and [8] we show the p/p ratio and p flux for four hadronic interaction mod- 
els and three CR propagation models. To get these result, besides the scale factor cp we 
adopt a charge-dependent solar modulation potential to fit the AMS-02 data better. The 
colored bands of the p/p ratio and the grey bands of p flux represent the uncertaities of 
the corresponding propagation models with the propagation parameters varying within 9596 
confidence ranges. Note that the GALPROP default hadronic interaction model underesti- 
mates antiprotons at low energies compared with the AMS-02 data for the DR propagation 
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FIG. 8. The same as [6] but for the DC2 model. 


model. 


a factor 7 in the DR-2 model [29]. Here we note that this underproduction actually depends 
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Such an underproduction has been realized for a long time when studying the 
PAMELA n" and BESS results bi . In order to solve this problem, the dif- 


fusion coefficient at low energies should be modified |77]. This is a motivation for introducing 


on the hadronic interaction models too. On the contrary, the QGSJET-II04m model would 
overproduce antiprotons at low energies in the DR propagation model. 

We find that the EPOS 1.99 model almost gives the best fit to the AMS-02 data in all 
three propagation models. It naturally predicts a flat p/p ratio at the high energy end. On 
the other hand, the p flux predicted by the QGSJET-II 04m model is not quite consistent 
with the AMS-02 data. The EPOS LHC model, which nearly best fit the accelerator data, 
seems to underestimate the p/p ratio at the high energy end with obvious discrepancies 
compared to last several data points. 

In Table [IT| we show the parameters that can give best fit of the p/p ratio to the AMS- 
02 data for the different propagation and hadronic interaction models. The propagation 
parameters are within the 95% C.L. range fitting to B/C. Here ó, is the solar modulation 
potential for protons by fitting to the AMS-02 proton spectrum. The solar modulation 
potential value for antiproton 6; is only slightly different from ¢,. We notice that the rescale 
factor cp is near 1 within 20%. It can be seen that the EPOS 1.99 and the GALPROP default 
models can well fit the AMS-02 data with small X? values. EPOS LHC also provides an 
acceptable fit to the data, while QGSJET II-04m tends to give a large x°. 


V. IMPLICATIONS FOR DARK MATTER ANNIHILATION 


As shown in the previous section, the theoretical prediction of the antiproton flux has 
quite large uncertainties. In some cases the prediction fit the AMS-02 data very well, which 
can give a strong constraint on the dark matter annihilation. However, in some cases, there 
are discrepancies between the prediction and the data, which may indicate a contribution 
from dark matter. We discuss the implications on DM annihilation according to the AMS-02 
antiproton result in this section. 

We use the results of PPPC 4 DM ID including the electroweak corrections to 
calculate DM annihilation. The DM density distribution is assumed to be the Navarro- 
Frenk-White (NFW) profile with the local density of 0.3 GeV cm~’. In Figs. [9] [10] 
and [II] we set upper bounds on the DM annihilation cross sections in the bb and WTW- 
channels at 9596 C.L. by using the AMS-02 p/p result. Note that if the discrepancies between 
the predicted p/p ratio and the AMS-02 data at the high energy end are taken seriously, 


there may be room for DM signals at high energies. Therefore, we also attempt to introduce 
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UA 


dV./dz 


Tan & Ng + BPO1 EPOS LHC EPOS 1.99 QGSJETII-04m 


DR DR? DC DR DR2 DC DR DR? DC DR DR2 DC 

115 1.05 0.99 0.899 0.959 0.76 0.923 0.853 0.836 122 12 104 

0.429 0.69 0.773 0.468 0.928 0.809 0.462 0.592 0.823 0.58 0.843 0.823 
77.9/57 38.9/57 51.6/57 131/57 78.1/57 300/57 119/57 28.9/57 57.3/57 119/57 298/57 1419/57 
0.515 0.745 0.644 0.562 0.773 0.674 0.554 0.71 0.686 0.483 0.703 0.686 
7.76 287 14 4.668 2.79 239 605 291 18 719 434 18 

0.337 0.433 0/0.453 0.338 0.42 0/0.478 0.335 0.418 0/0.523 0.324 0.39 0/0.523 
= = Aa = 302 = å x Bre = = Fe 

546 28 169 492 2.78 313 438 286 217 479 392 217 

41.3 2909 — 344 189 — 33.9 195 — 425 244  — 

= = 0826 = = 29 = = 268 — — 268 


a The ô would be 0 for R < Ro in the DC case 


TABLE II. The values of parameters for the best fit to data in different scenarios. 


a DM signal to improve the fitting quality, and show corresponding regions of the parameter 


space in Figs. D] and [I] 


Fig. S]shows the upper limits on the DM annihilation cross sections for the bb and WtW- 
final states for the QGSJET II-04m model. Here we only give the constraint in the DR 


propagation model, since in the other two propagation models the secondary p/p predictions 


can not fit the data to derive a reasonable bound. The red bands represent the uncertainties 


from the propagation model. For comparison, the limits from the Fermi-LAT gamma-ray 


observation of dwarf galaxies are also shown in the plots. We can see that the constraints 


from the AMS-02 p/p result are more stringent than the Fermi-LAT constraints. We also 


find a small parameter region located around 200 GeV, where a DM annihilation signal 


would improve the fitting quality. 


In Fig. [10] we set constraints for the EPOS-LHC model. The constraints are worse than 


the Fermi dwarf galaxy constraints, since the EPOS LHC model generally predicts a slightly 
lower p/p ratio than the data at the high energy end. Recalling that the EPOS LHC 


model fits the accelerator data well, such discrepancies should be taken seriously. In order 
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FIG. 9. The 95% upper limits on the dark matter annihilation rate (ov) derived from the AMS-02 


p/p ratio for QGSJET II-04m with the DR propagation model. The left panel is for the bb channel, 


and the right panel is for the W*W~ channel. The red bands indicate the propagation uncertainty 


in 9596 C.L., while the red line is corresponding to the case that fit to the p/p best. For comparison, 


we also show the upper limits from the Fermi-LAT observation of dwarf galaxies " as the blue 


lines. The contour denotes a DM signal in the 6896 confidence region favored by the AMS-02 p/p 


ratio data in the m, — (cv) plane. 
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FIG. 11. The same as Fig. D| but for EPOS 1.99 interaction model with the DR-2 and DC 


propagation models. 


to improve the fit, a DM signal with a mass of ^ TeV can be introduced. However, the 
corresponding DM annihilation signal seems to be excluded by Fermi. Here we mention that 
such a contradiction may be evaded, if a velocity-dependent DM annihilation cross section 
is introduced . 

We give the constraints for the EPOS 1.99 model in Fig. As the EPOS 1.99 model fits 
the AMS-02 antiproton data best, it sets very stringent constraints on the DM annihilation 


cross section. 


VI. CONCLUSIONS 


In this work we recalculate the astrophysical CR p/p ratio and p flux, and compare 
them with the latest AMS-02 results. We investigate the uncertainties from the propagation 
models and the hadronic interaction models. 

We consider the DR, DR-2, and DC propagation models, where the prorogation param- 
eters are determined by fitting the B/C and !?Be/?Be data with a MCMC algorithm. All 
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these models can provide a very good fit to the B/C data. 

We note that the uncertainties from propagation and hadronic interaction may interplay 
with each other. For example, it has been recognized that the DR propagation model may 
underestimate the p flux at low energies by GALPROP bz 78 . However, we find that 
the p flux is even overproduced at low energies if the QGSJET II-04 model is adopted. For 
the DR-2 and DC propagation models, the predictions given by QGSJET II-04m seem not 
consistent with the AMS-02 p data, while other three interaction models can fit the data 
quite well. 

We further derive the upper bound on the DM annihilation cross section from the AMS- 
02 p/p data. It is found that in most cases the bound can be stronger than the constraint 
from Fermi observation of the dwarf galaxies. 

Among these interaction models, the EPOS LHC model can fit the NA49 and LHC data 
very well. However, this model predicts a slightly lower p/p ratio than the AMS-02 data 
at the high energy end. If such a discrepancy is taken seriously, a ~ TeV DM annihilation 
signal is necessary to fit the AMS-02 data. 
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